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Abstract 

The Regge-Mueller formalism is used to describe the inclusive spectrum of the 
proton in pp collisions. From such a description the energy dependences of both 
average inelasticity and leading proton multiplicity are calculated. These quantities 
are then used to establish the connection between the average charged particle 
multiplicities measured in e + e~ and pp/pp processes. The description obtained for 
the leading proton cross section implies that Feynman scaling is strongly violated 
only at the extreme values of xp, that is at the central region (xf ~ 0) and at the 
diffraction region {xp ~ 1), while it is approximately observed in the intermediate 
region of the spectrum. 
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1 Introduction 



It is experimentally well known that the energy dependence of the charged particle mul- 
tiplicities in e + e~ and pp/pp processes exhibit a quite similar behavior. In the late 70's, 
experiments analysing pp collisions at the CERN ISR Collider have shown that not 
the total center-of-mass energy y/s is used for particle production; instead, a considerable 
fraction of the available energy is carried away by the leading proton. These experiments 
have shown that a more adequate way of comparing average multiplicities from different 
reactions is in terms of the amount of energy effectively used for multiparticle production. 
The problem is how to determine this quantity. 

Observations like these have inspired several attempts to describe (n ch ) e+e _ and (n c h) pp 
by an universal function. In ref. @, for instance, two corrections are made to compare 
these quantities: the energy variable for (n c h) pp is corrected by removing the portion 
referring to the elasticity (the fraction of the energy taken by the leading particle) and 
then the average leading proton multiplicity is subtracted. A similar idea is followed in 
ref. || where attempts are made to establish this universal behavior by fitting. 

In the present paper, we analyze the same subject by following an analagous point of 
view, but rephrasing the procedure in the following way. It is assumed that, if in e + e~ 
collisions the average charged particle multiplicity is given by 

{n ch ) e+e . = tf(>/5), (1) 

then in pp/pp collisons we have 

(n ch ) pp = (n ) + N((k p )^), (2) 

where N(W) is an universal function of the energy available for multiparticle production, 
W, (n ) is the average leading particle multiplicity, and (k p ) is the average inelasticity. 
In and 0], the quantities related to (no) and (k p ) are supposed to be constant. In 
particular, in || they are determined by a simultaneous fit of (n ch ) e+e _ and (n ch ) pp data. 

Our procedure, instead, consists in obtainning these quantities ((no) and (k p )) not 
from fitting (n c h) e+e - and (n c h) pp , but in a totally independent way, from the inclusive 
reaction pp — > pX, paying particular attention to their energy dependence. After doing 
that, the obtained (no) and (k p ) are applied to (0) via a parametrization of ([[]) and the 
result is compared to data in order to verify to what extent such a hypothesis is acceptable. 

This procedure seems to be very well defined and straightforward, but it should be 
noticed that it drives to some difficult problems. The question is that it requires a previous 
knowledge about the energy dependence of the inelasticity and about the behavior of the 
average leading particle multiplicity which constitute themselves problematic subjects. 
In particular, the energy dependence of the average inelasticity is a very much disputed 
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question since there are opposite claims that this quantity increases ||, ||, || [F] or that 
it decreases || || [K]] with increasing energy at quite different rates. In spite of the of 
models predicting extreme behaviors, i.e. very fast increase of the inelasticity (like in 0) 
or very fast decrease (like in ||), most of these analyses referred here report the average 
inelasticity as having a smooth and slowly changing behavior. [] This is once again verified 
here in a different way. 

The idea of discussing the energy behavior of the average inelasticity in connection 
with the energy dependence of (no) and (k p ) is not new. Of particular interest to present 
work is an analysis with this purpose performed by He |10[ . He has extracted values of the 



average inelasticity by using arguments similar to those given above and obtained results 
pretty much in agreement with the predictions of ref. ||. We shall argue below that such 
an agreement is probably due to the fact that two important effects are missing in his 
analysis. 

Another controversial question involved in the present analysis (but treated here just 
en passant) is that referred to unitarity violation in diffractive dissociation processes. 
This is an old-standing problem that has come back to the scene due to the fact that 
recent measurements on hard diffractive production of jets and W's revealed a large 



discrepancy between data and theoretical predictions. In ref. []12f, it is proposed that such 
a discrepancy in hard diffraction has to do with unitarity violation in single diffractive 
processes. Since we are going to describe the inclusive reaction pp — > pX, we have to face 
this problem in the region of the spectrum where diffractive processes are dominant. 

A by-product of the present analysis is a complete parametrization for the reaction 
pp — > pX in the whole phase space. This is obtained basically within the Regge- Mueller 



approach [13], but including the modifications suggested in [12] for the diffractive contri- 
bution. 

The paper is organized as follows. In Section 2 we present the theoretical framework 
used to describe the leading particle sprectrum. Section 3 is devoted to show how this 
formalism is applied to describe the experimental data. In Section 4 we discuss the con- 
nection between (n ch ) e+e _ and (n ch ) pp . Our main conclusions are summarized in Section 
5. 



2 Leading particle spectrum 

In order to perform our analysis, we need to calculate the quantities 

(n ) = [ ~r~~ dx F (3) 

&inel J d^F 



x For a recent account on this subject from the viewpoint of cosmic-ray data, see ref. [p"T| . 
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and 

(A; p ) = 1 — (x.f) = 1 / x F - — dx F (4) 

&incl J dxp 



as a function of energy. We apply the Landshoff parametrization [FJ] <Ji ne i{s) = 56 s + 
18.16 s [mb] to represent the inelastic cross section within the energy range where 
multiplicity data are included, and the inclusive cross section is simply given by Oi nc i = 
J dxp. Thus, the whole analysis depend on the knowledge of the leading particle 
spectrum dafdxp and its evolution with energy. The obtainment of this spectrum is 
detailed in the discussion that follows. 

The invariant cross section for the inclusive reaction ab — > cX is given by 

d 3 a 2E d 3 a 

where xp = Ippj^fs is the Feynman variable for the produced particle c and E, pp, pp 
are respectively its energy, longitudinal and tranversal momenta. Particularly in the 
diffractive region (xp ~ 1) such a quantity is usually expressed in terms of 

£,d 3 cr s d 2 a Xp d 2 a 
dp 3 7r dt dM 2 tcx f dt ' 

with xp = 2E/y/s, £ = M 2 /s = 1 — xf and —t = m 2 {1 — xf) 2 fxp + p\jxp. Variable 
M 2 is the missing mass squared defined as M 2 = (p a + pb — p c ) 2 - 

The procedure to calculate the invariant cross section employed here comes from the 
Regge-Mueller formalism which consists basically of the application of the Regge theory for 
hadron interactions to the Mueller's generalized optical theorem. This theorem establishes 
that the inclusive reaction ab — > cX is connected to the elastic three-body amplitude 
A(abc — > abc) via 

E—^-(ab — > cX) ~ -Disc M 2 A(abc — ► abc), (7) 
ap^ s 

where the discontinuity is taken across the M 2 cut of the elastic amplitude. It is assumed 
that this amplitude in turn is given by the Regge pole approach. Different kinematical 
limits imply in specific formulations for the invariant cross section at the fragmentation 
and central regions. In the following, we specify the concrete expressions that these 
formulations assume in such regions (details can be found in fD| ) . 

A. Fragmentation Region 
In our description, the invariant cross section for the reaction pp — > pX at the frag- 
mentation region is compounded of three predominant contributions which are determined 
within the Triple Reggeon Model (this is the particular formulation that (0) assumes in 
the beam fragmentation region with the limits M 2 — > oo and s/M 2 -> oo |T§). These 
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contributions, depicted in Fig. 2, correspond to pomeron, pion and reggeon exchanges and 
are referred to as IP IP IP, 7T7rIP, IRIRIP, respectively. 

In the diffractive region, the IPPIP contribution is dominant and (we assume for the 
reasons given below) is given by 



' d 2 a s 
dtd^ 



IPJPIP 



where fw,Ren(£,,t) is the renormalized pomeron flux factor proposed in [0 with the pa- 



rameters defined in JT5|, that is 



/p,Ren(£>*) - 
with the Donnachie-Landshoff flux factor [16| 



N(s) 



and 



N(s) 



IQtt 



i rt=o 



f F {£,t) dtdt 



'l.5/s J — oo 

In the above expressions, F\(t) is the Dirac form factor, 



Fi(t) 



1 



(4m 2 - 2.79*) 
(4m 2 -t) (l-^r 



(9) 



(10) 



111 



(12) 



the pomeron trajectory is aijp(t) = 1 + e + at with e = 0.104, a = 0.25 GeV~ 2 and 
/3q = 6.56 GeV -1 , determined from |T7[. In Eq.(|8]), the pomeron-proton cross section is 
given by 

a Pp (M 2 ) = ft gw (sO e (13) 

with the triple pomeron coupling determined from data as g-p = 1.21 GeV^ 1 . 

Since this scheme to calculate the diffractive contribution is not the usual one, some 
comments are in order. The usual derivation of the Triple Pomeron Model gives ([10|), 
the standard flux factor, instead of (|]), the renormalized one. The problem is that the 
standard flux factor drives to strong unitarity violation and the renormalization procedure 



was conceived [|T^] as an ad hoc way to overcome this difficulty. Although a rigorous 
demonstration of the renormalized scheme is still missing, it is acceptable in the sense 
that it provides a good description for the experimental data at the diffractive region (see 
a detailed discussion in fl8fl ). 

The pion contribution (tttiJP) is given by |L9 
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dtd£ 



ttttF 



where 



47T 47T (t — p : 



2\2 



^(*-p 2 )^l-2a vr (t) 



(14) 



(15) 



and a n (t) = 0.9(t - p 2 ) with p 2 

constant in g 2 /47r = 15.0 and putting &„- = (see also 
<t*p(sO = 10.83 (sO - 104 + 27.13 (sO -0 - 32 [m&] is taken from fL7 



0.02 GeV . We follow [EtJ in fixing the coupling 
The pion-proton cross section 



If one considers only the diffractive and near-to-diffractive regions and low px (— t ~ 
0.0 — 0.1 GeV 2 ), the contributions outlined above are enough to provide a good description 
of the available data 



sec 
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However, when one wants to consider larger and 
xp < 0.9, at least a third contribution is required. That is the reason why we introduce 
the reggeon contribution. 



The reggeon contribution (1R1RIP) is determined by 

/ d 2 a\ 



x 0-irpOO 



IFURIP 



with 



and 



fn(€,t) = -r^e 

l07T 



(16) 



(17) 



o"]Rp(sO = An ^n(sO e - (18) 

In this case, the trajectory is assumed to be an(t) — 0.5 + 1 while the constants (3-^ = 
(Pom. 9m) an d 6r remain to be determined from data. 

Thus, with the expressions and parameters given above, the IPPIP and tttcJP contri- 
butions are completely specified; only the 1R1RIP contribution remains to have the final 
parameters determined. 

B. Central Region 

In order to describe the leading particle spectrum in the central region, we use the 
Double Reggeon Model [O] that gives the invariant cross section as 



E 



d s 



a 



dp 3 



E 7y' 



hi 



t 


Oi(0)-1 


U 


Qy(0)-1 


so 




so 





(19) 



where my = {p\ + m 2 ) 1 / 2 is the transversal mass, and u = —rriT^/s e y and t = 
—m<r\fs e y are the Mandelstam variables given in terms of the rapidity y = In- 1 ' 



rriT 
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Function ^fij{mj) corresponds to the product of the three vertices of the diagrams de- 
picted in Fig. 3. These diagrams represent the contributions taken into account in the 
present analysis: IP IP, IPIR + IRJP, and IRIR (pion contributions are not considered in this 
case because they are totally covered by the others). 

We assume for the coupling function 7r,(my) a simple gaussian form, 

ll3 {m 2 T ) = Yij e~ a - m T (20) 

where is a constant that already embodies the product of the couplings belonging to 
the triple and quartic vertices. With these definitions, the invariant cross sections for the 
three contributions become 





= 2 r P]R e - a ™ m T ( TOT ^y+^m-i cosh ^ + e _ Q]R (o)) y ], (22) 

P1R+EP 



and 




= r R ia e~ a ^ m T (m T ^^("l- 1 '. (23) 

KIR 

In the above expressions again aja(0) = 0.5 and e = 0.104 |17|]. Differently from the 
fragmentation region where almost all parameters are already established, in this region 
almost all of them (expect for the intercepts just mentioned) must be fixed from data. 

The expressions given above could be enriched by detailing the reggeon exchange in 
terms of /, p, u, a 2 , and taking into account all crossed terms, but in fact we are pursuing 
here a minimal description in which only the dominant and effective contributions are 
considered. We shall see below that these contributions are enough to provide a good 
description of the available data. 



3 Description of experimental data 

Experimental data on leading particle spectrum are very scarce. A compilation for pp — > 
pX is shown in Fig.l where data from three experiments |T|, |2T], |35| are put together (the 
curve and the insert in this figure should be ignored for the moment). As can be seen, 
a pretty flat spectrum is exhibited, except for xp ~ 1 where the typical diffractive peak 
appears.^ 

2 This peak is absent from the Aguilar-Benitez et al. data due to trigger inefficiency for xp > 0.75 in 
this particular experiment pf| . 
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The problem that arises when one tries to describe the pp — > pX reaction in the whole 
phase space is that the available data are not enough to determine unambigously each 
one of the contributions outlined above. One may have noted in the previous section 
that we have summarized all secondary reggeon exchanges (except for the pion) in a 
single contribution denoted by 1R and the reason is the following. When one analyzes, 



for instance, total cross section data (like in |P7|), it is possible to establish (to a certain 



extent) the relative amount of the different contributions. Actually, this is enforced by 
the changing shape exhibited by the data in different regions. That is not the case here 
because out of the diffractive region the spectrum is pretty flat and that makes it difficult 
to discriminate the regions where the different exchange processes contribute the most. 
Thus, in order to establish how the expressions outlined above are summed up to compose 
the observed spectrum, we have to follow a particular strategy. 

Since our intention was obtainning an acceptable description for pp — > pX data in the 
whole phase space, we did not use in our fitting procedure the data shown in Fig.l which 
represent only the x^-dependence. Instead, we have set those data apart to be used only 
at the end to check our final results which, in fact, were obtained with distributions giving 
in terms of both xf and pr dependences. 

Our procedures to determine the contributions at the central and at the fragmentation 
regions are quite different. The main problem is that these regions overlap each other and 
thus it is pratically impossible to separate them (or establish clear limits). To overcome 
this difficulty we assumed that, except for normalization effects, the xp and pt depen- 
dences of the proton produced in the central region through the reaction pp — > pX is the 
same as for the antiproton produced in pp — > pX. This assumption was implemented by 




central 



fitting simultaneously the data shown in Figs. 4 and 5 |23|, 23] through the expressions 

(24) 

\ a f / pp->pX \ / PP \ «P / PE+RP \ «P / Kir 

and 

central / ,q \ central 






(25) 

pp—>pX \ "** / pp— >pX 

The idea is that the data of Fig.4 provide the information on the xf and px dependences 
through Eqs. (|2l"|)-(|24|) and relation xf = 2rriTsinh(y) j \/s, while the connection between 
pp — > pX and pp — > pX is established by fitting the data of Fig. 5 through the function 
X(s) of Eq.(|25|). The parameters and a^- of this fit are given in Table [I] while A(s) is 
parametrized as A(s) = 1.0 + 11.0 s~°' 3 . 

The agreement with data of Figs. 4 and 5 is not perfect, but that is because we are 
simplifying the description by considering only a few contibutions, the dominant ones. As 
stated before, this is enough for the purposes of the present analysis. 
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Now we are able to obtain the total description by adding up central and fragmentation 
region contributions. As explained before, the contributions dominant at the fragmenta- 
tion region, Eqs. (PD-flTSP, are almost completely determined. The parameters /3jr and 
6]r referring to the 1R1RIP contribution are established by fitting the data of Fig. 6 (from 



22|). This is done by using the expression 



t q \ total / 7 q \ / t q \ / ,q \ / 7 q \ central 





^P 3 /ppp V dp 3 J nnJP V ^pVfjrp \ dp 3 ) pp _ >pX 

where the last term refers to Eq. (^) with the parameters given in Table [l|. With this final 
fit the remaining parameters result to be (3jn = 2465.7 mb GeV~ 2 and = 0.1 GeV~ 2 . 
Fig. 7 offers a view of how the different contributions are composed to form the final result 
and how this picture evolves with p?. 

The different contributions of the invariant cross section in both regions integrated 
over pt produce the results of da/dxF for both reactions exhibited in Fig.l (solid curves) 
for pi ab = 400 GeV/c. We remind the reader that these data were not used in the fit, 
but are used now to check the reliability of the whole procedure. From this figure it is 
possible to see that the final description obtained for the leading proton spectrum is quite 
reasonable. 



4 Connection between (n c h) e + e - and (n c h) pp 

The results obtained above specify completely the behavior of the leading particle spec- 
trum and allow us to calculate (n ) and (k p ) as given by (|) and ([|). In Fig. 8, we show the 
energy dependence of these quantities as obtained in the present analysis (solid curves). 
In the same figure, it is also shown the average inelasticity as predicted by the Interact- 
ing Gluon Model (IGM) || (dot-dashed curve) for comparison. The average inelasticity 
obtained from the present analysis is very slowly increasing with energy, close to the 
behavior predicted by the Minijet Model 0. 

With these results we can come back to our original intent which is checking the 
hypothesis of universal behavior of the multiplicity that is specified by Eqs.([L|) and (|2|). 
In order to do that, we first establish a parametrization for N(y/s) through 

N(y/s) = at + a 2 ln{—) + a 3 ln 2 (-) (27) 

with so = 1 GeV 2 . However, before performing the fit to experimental data, an additional 
effect has to be considered. This is because, besides the charged particles produced at 
the primary vertex, (n c h) e+e - data include also decay products of K® — > 7i + n~ , A — > pn~ , 
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and A — > pir + . Following ||, we take this contamination into account by computing the 
ratio R = (n c h) k °,a,a/ ( n ch) e+e - and by redefining (P as 

M e+e - = y^, (28) 

with R = 0.097 ± 0.003. This value was taken from || and, besides the references 
quoted therein, it is in good agreement with experimental data from ref. No energy 
dependence for R can be inferred from these data. The fit using (|27| ) and (^) gives 
ax = 2.392, a 2 = 0.024, and a 3 = 0.193. 

In Fig. 9, we show the above parametrization describing (n c h) e + e - data from references 
quoted in |26) and the calculated curve for (n c h) pp in comparison with data from The 



agreement with these data enables us to consider that our premises about the universal 
behavior of (n c / l ) e+e „ and (n c h) pp are confirmed. Of course, this conclusion is restricted 
to the energy dependence of (n ) and (k p ) shown in Fig.8. 

The solid curve of the insert in Fig. 9 shows what happens when the IGM average 
inelasticity is applied to the same purposes. One could argue that this last result is 
conditioned by the use of (no) obtained in the present analysis which increases with 
energy. However, we note that increase in (n Q ) plays against increase in (k p ) since these 
are competitive effects. 

Now a comment on the He analysis |H|, where the relation 



nt h e (y/s e + e -) = n^(k(y/s^)y/s^) (29) 



b ch 

is employed. After fitting and n v ^ h independently, He imposes that relation ( p9"D 

holds and extracts the inelasticity k from this assumption. This is similar to what we 
have done, but we think that the result of decreasing inelasticity and the agreement with 
IGM obtained in such an analysis comes from the fact that neither the leading particle 
multiplicity (no) nor the effect of decay products (R) is considered and we see no reason 
for ignoring such effects. 

A surprising outcome of the present analysis is shown in Fig. 10 (a) where the normal- 
ized cross section l/<Ti nc i da/dxF is calculated up to the LHC energy. It is shown that, if 
the present description holds up to such high energies, Feynman scaling is approximately 
observed in the intermediate fragmentation region, 0.2 < xf < 0.8, but is violated in op- 
posite ways at the central and diffractive regions. Fig. 10 (b) shows the same results but 
in a scale that makes more evident the scaling violation at the central region. This result 
seems to say that the increase of production activity at the central region occurs at the 
expenses of a supression of the diffractive processes. However, this is just a speculative 
observation that should be investigated more thoroughly. 
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5 Conclusions 



We have presented in this paper a description of the inclusive reaction pp — > pX in the 
whole phase space within the Regge- Mueller formalism, modified by the renormalization 
of the diffractive cross section. The average multiplicity and the average inelasticity were 
obtained from the leading proton spectrum and both of them resulted to be increasing 
functions of energy, in agreement with [B, |5| and particularly with J7|. The energy depen- 
dence of these quantities is such that allows one to accommodate very well the charged 
particle multiplicities (n ch ) e+e _ and (n ch ) pp by an universal function once an appropriate 
relation is used. 

An additional result is that the normalized leading proton spectrum approximately 
observe Feynman scaling for intermediate Xp, whereas such scaling is violated at the 
central and diffractive regions. 
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u 


Tij (mb GeV~ 2 ) 


(Hi (GeV- 2 ) 


PIP 


23.53 


3.90 


P1R 


-29.8 


3.45 


KJR 


13.75 


1.80 



Table 1: Values of the parameters and a^. 
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Figure 1: Inclusive spectrum for the reactions pp — > pX and pp — > pX (in the insert). 
Data from [[|, [H], The solid curves are the results of the fit described in the text 
calculated for 400 GeV/c. 
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Figure 2: Triple-reggeon diagrams considered in the present analysis for the reaction 
pp — > pX. The particles corresponding to the external lines are all as in the first diagram. 
These diagrams represent the contributions that are dominant in the fragmentation region. 
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Figure 3: Double-reggeon diagrams considered in the present analysis for the reaction 
pp — > pX. The particles corresponding to the external lines are all as in the first diagram. 
These diagrams represent the contributions that are dominant in the central region. 
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Figure 4: Invariant cross section for the reaction pp — > pX at the ISR energies. The 
description is obtained with Eqs.(f21|)-(|24|) and parameters of Table 1. Data taken from 
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Figure 5: Invariant cross section for the reaction pp 
obtained with Eq.fl25|). Data taken from |24 . 



pX at xf = 0. The description is 
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Figure 6: Invariant cross section for the pp — > at the fragmentation region. Curves 
calculated with Eq. (|26|) . Data from |22 . 
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Figure 7: Invariant cross section and its contributions for pp — > pX at y^s = 53 GeV 
for two pr values. These plots show how the interplay among the different contributions 
changes as pr increases. 
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Figure 8: Energy dependence of (a) average leading proton multiplicity (no) and (b) 
average inelasticity (k p ). In the lower figure (b) it is shown (k p ) obtained in the present 
analysis (solid curve) compared to the same quantity as predicted by the IGM || (dot- 
dashed curve). 



19 



50 



40 



30 



A 

O 

V 



20 



10 



~~ i 1 r~ 



- 1 1 — i — r - 



10" 









/ / ■ 








M / 

w / 




M ] 



10' 10 




Oee 
App 
■ pp 

fitted curve 

Calculated curve 



i i i i 



i i i i i i i 1 1 



i i i i_ 



10' 



io 2 

Vs (GeV) 



10 J 



Figure 9: Main figure: average charged particle multiplicities (n c / l ) e+e _ and {n c h) pp as a 
function of the center-of-mass energy. The dot-dashed curve refers to the fit obtained with 
Eqs. (p7[) and ([281). The solid curve was calculated with Eq.(^) by using (no) and (k p ) as 
calculated in the present analysis. Insert: the same as in the main figure, but using the 
average inelasticity given by the IGM M. 



20 





i\ i , , , , — 1 1 — 

u -1 -2 -1 

10 10 10 10 

x F 

Figure 10: (a) Normalized leading proton spectrum calculated up to the LHC energy, (b) 
The same as (a) but with logarithmic scale for x F . 
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